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1.0  PUBLICATIONS 


During  the  year  a  number  of  publications  appeared 
whose  research  was  sponsored  by  the  Grant.  Copies  of  the 
first  page  of  the  papers  are  attached  and  will  indicate  the 
Journals  in  which  the  papers  appeared  and  the  abstracts  will 
give  an  indication  of  the  research  contained  in  the  papers. 


2.0  RESEARCH  DIRECTIONS 


In  addition  to  the  research  reported  in  the  pub¬ 
lished  papers  referred  to  above,  some  progress  was  made  in 
relating  the  inverse  scattering  problem  to  causality.  In 
particular,  Newton’s  "miracle”  formula  for  the  potential  in 
three-dimensional  inverse  scattering  was  derived  purely  from 
causality  considerations.  Moreover,  the  one -dimensional 
analogue  indicated  some  possible  problems  when  point  eigen¬ 
values  were  present.  Termination  of  the  grant  precluded 
completion  of  the  paper  and  subsequent  publication. 
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Eigenvalues  and  eigenfunctions  associated  with  the  GePf and-Levltan 
equation 

Harry  E.  Moses*' 

Center  for  Atmospheric  Research,  University  of  Lowell  Lowell,  Massachusetts  01854  I 

Reese  T.  Prosser  j 

Department  of  Mathematics.  Dartmouth  College.  Hanover,  New  Hampshire  03 753  I 

(Received  8  June  1983;  accepted  for  publication  10  August  1983)  j 

It  is  shown  here  that  the  solutions  of  the  Gel'fand-Levitan  equation  for  inverse  potential 
scattering  on  the  line  may  be  expressed  in  terms  of  the  eigenvalues  and  eigenfunctions  of  certain 
associated  operators  of  trace  class.  The  details  are  sketched  for  the  case  of  rational  reflection 
coefficients,  and  carried  out  for  the  simplest  class  of  examples. 

PACS  numbers:  03.80.  +  r,  03.65.Nk 


1.  INTRODUCTION 

The  Gel'fand-Levitan  equation  plays  a  central  role  in 
solving  inverse  scattering  problems  in  one  dimension. 1  In  the 
case  where  the  problem  involves  a  scattering  potential  K(x) 
defined  for  —  oo  <x  <  ■+  oo ,  for  example,  we  know  that 
V{x)  may  be  recovered  from  the  reflection  coefficient  r\k ), 
defined  for  —  oo  <k<  +  oo,  as  follows:  set 


it  follows  that  the  operator  P  (w)R  is  of  trace  class  for  each  w, 
and 

j 

tr  f»(u>)R«J*  ?(2 i)dz.  (9) 

One  can  then  define  the  Fredholtri  determinant  A  (u>)  of  the 
operator  (/  +  P(w)R )  by  (cf.  Ref.  :3,  p.  2351!) 


R(x,y)  =  ?(*+y)  =  -^J  + V^**"^**.  (1) 

and  then  solve  for  K  (x,  y)  the  Gel’fand-Levitan  equation 

R>,.y)  +  /t(x,y)  +  j  K{xj)R\?,y\dx=*0.  (2) 

Then  the  potential  V{x)  appears  as 

K(x)  =  2 -£-*(*.*).  (3) 

dx 

(See  Ref.  2  for  a  general  discussion  of  this  procedure.) 

In  order  to  study  the  behavior  of  the  solutions  of  (2),  it  is 
useful  to  consider  the  associated  equation,  to  be  solved  for 
K{x,y,w): 

R  (x.y.w)  +  R  (x,y)  -f  J""  K{xj,w)R  {z,y)di~0.  (4) 

Evidently  K {x, yjc)  =  K\ x, y).  Now  (4)  may  be  expressed  in 
operator  form  with  w  as  a  parameter: 

K(w)  +  R  +  K{w)P(w)R  »  0.  (5) 

Here  R,  K (w),  and  P(w)  are  integral  operators  with  kernels 
R  (x,  y),  K  {x,  y.w),  and  P[x,  y.w),  with 

Plx,y,w)  =  0{w-x)6{x-y).  (6) 

Here  0(z)  is  the  Heaviside  function,  and  6{s)  its  derivative. 
Now  (4)  yields 


A  (tv)  ■>  det (/  +  P(w)JI ) 


«  exp  tr  log(7  +  J*(iv)R ). 

(10) 

Evidently 

log  4  (10) » tr  log(/  4-  P(w)R ) 

(ID 

and  so  . 

dw  A  (u>) 

- \x  P'(w\R[f  +  P\w)R )“• 

as  —  \tP’(w)R\LW)- 

(12) 

Here  we  have  used  (8).  But  P'(w)R  (ui)  has  kernel 
6{w-x)Klx,yjv),ao 


-  trP*(i»MC(u>)  *  -  J"  6{w  -  x\K  \?jc,w\dx 
=  —K  \w,w,w) 

-  —  K  (u;,u»).  (13) 

Hence  by  (3) 

V(w)~2~K[w,w) 

dw 

«- 2-^1  log  J(«*  (!♦) 


R  (w)(J  +  P[w)R  )=—/?,  (7) 

and  hence,  whenever  (/  +  P(w)R )  is  invertible, 

A(u/)  =  -R{I+P(w)R)-'.  (8) 

Now  suppose  that  the  reflection  coefficient  t{k )  is  such 
that  its  Fourier  transform  ?(z)  is  smooth  and  integrate.  Then 
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This  formula,  which  gives  V  directly  in  terms  of  R,  first 
appears  in  Ref.  4,  and  has  since  been  rediscovered  by  several 
authors,  including  us.5  In  one  sense,  this  formula  by-passes 
the  Gel'fand-Levitan  equation,  since  it  gives  V  directly  in 
terms  of  R,  and  once  V  is  known  everything  about  the  scat¬ 
tering  problem  is  known,  at  least  in  principle. 

In  another  sense  ( 14)  is  no  better  than  (4),  since  the  cal¬ 
culation  of  the  determinant  d  (u»)of(/  +  P(w\R )  is  not  usually 
an  easy  matter  in  practice.  One  possible  approach  is  to  calcu- 
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For  the  moat  part,  the  phase-retrieval  problem  has  been  dominated  by  the  use  of  the  Hilbert  transform  on  the  loga¬ 
rithm  of  the  absolute  value  of  the  wave  function.  In  our  approach,  by  contrast,  we  use  the  intensity  associated  with 
the  function  and  the  intensities  in  Fourier-  or  Mellin-tramform  space  corresponding  to  a  set  of  ape:  lures  to  obtain 
the  phaaes  within  a  constant  by  means  of  simple  formulas.  The  set  of  apertures  required  is  those  arbitrarily  close 
to  a  fired  aperture  and  its  complement.  The  higher-dimensional  cases  are  also  treated  for  the  case  of  the  Fourier 
transform. 


fffTRODUCnON 

TW  Boat  commonly  used  method  used  to  obtain  the  phase 
•fa  complex  function  f(x )  of  a  real  variable  x  it  to  take  the 
fagarithm  of  /(x).  which  in  certain  circumstances  can  be  an¬ 
alytically  continued  in  the  complex  plane.  The  Hilbert 
transform  is  then  used  to  find  the  phase,  which  is  the  ratio  of 
the  imaginary  part  of  the  logarithm  to  the  real  part;  the  real 
part  ■  the  logarithm  of  the  absolute  value  of  the  function.  For 
unique  results,  the  xeroa  of  the  analytic  function  must  be 
prescribed.  Various  methods  are  used  to  find  them  from 
physical  measurements,  which  correspond  to  the  use  of  ap¬ 
propriate  sets  of  apertures.  A  survey  of  approaches  to  the 
phase-retrieval  problem  is  given  in  Ref.  1.  A  recent  method 
hr  locating  the  reroe  is  given  in  Ref.  2. 

As  far  as  we  can  determine,  the  first  departure  from  the  use 
at  the  Hilbert  transform  was  given  in  1963  by  Lomont  and 
Moses.3  The  object  of  their  approach  was  to  give  data  in 
terms  of  intensities  only.  It  was  shown  that  a  necessary  and 
sufficient  condition  for  two  complex  functions  with  the  same 
amplitude  to  have  the  same  phase  within  a  constant  was  that 
the  squares  of  the  amplitudes  in  Fourier-tranaform  space  be 
equal  for  ail  apertures  through  which  f(x)m  observed.  In  1971 
Gerchberg  and  Saxton4  independently  used  the  input  required 
by  the  theorem  of  Ref.  3  to  obtain  the  phase  and  gave  nu¬ 
merical  methods  to  obtain  it.  In  Ref.  5  the  intensities  in 
Fourier-transform  space  are  also  used,  together  with  analytic 
properties  of  the  function  and  its  transform,  and,  like  the  work 
of  Ref.  4,  that  of  Ref.  5  goes  in  a  different  direction  from 
ours. 

In  this  paper  we  give  simple  formulas  for  phase  retrieval 
Used  on  measurements  using  apertures  arbitrarily  ekiae  to 
a  fixed  aperture  and  its  complement.  We  believe  that  our 
work  is  easier  to  use  in  many  cases  than  that  done  earlier,  al¬ 


though  numerical  implementation  has  yet  to  he  done  using 
our  method.  One  of  us  (Moses)  reviewed  the  material  of  Ref. 
3  from  a  different  point  of  view  and  showed  also  that  the  in¬ 
tensities  in  the  Mellin-transform  space  could  be  used  instead 
of  those  in  the  Fourier-tranaform  space.6 

Without  going  into  detail,  we  report  that  the  squares  of  the 
amplitudes  of  Fourier  transforms  correspond  to  intensities 
of  the  spectral  components  in  f(x)  when  the  signal  passes 
through  gratings  of  varying  aperture.  Another  well-known 
physical  situation  in  which  the  function  and  its  Fourier 
transform  occur  is  that  of  optical  systems,  in  which  it  can  be 
shown  that  the  wave  functions  in  the  aperture  and  the  image 
plane  constitute  a  Fourier-transform  pair.* 

The  applications  of  phase  retrieval  are  many  and  important. 
Among  the  applications  are  (dark-field)  electron  microscopy 
and  holography.  Other  applications  are  mentioned  in  Refs. 
1  and  5. 

The  significance  of  the  Mellin  transform  in  terms  of  scaling 
ia  discussed  in  Ref.  7,  in  which  a  method  is  also  given  whereby 
the  intensities  in  the  Mellin-transform  space  (in  one  dimen¬ 
sion)  can  be  obtained  as  a  Fourier  transform  of  a  scale  auto¬ 
correlation.  The  optical  importance  of  the  Mellin  transform 
is  discussed  in  Refs.  8  and  9,  in  which  methods  of  imple¬ 
menting  the  transform  optically  are  discussed. 

Although  the  Mellin  transform  can  be  obtained  from  the 
Fourier  transform  by  a  change  of  variables,  the  two  transforms 
represent  different  situations.  In  Ref.  7  it  is  shown  that  peaks 
in  the  Mellin-transform  power  spectrum  correspond  to  fea 
tures  in  the  original  function  that  are  periodic  in  mannifira- 
lion.  By  contrast,  peaks  in  the  Fourier-transform  power 
spectrum  correspond  to  features  that  are  periodic  in  (ran*- 
lat ion.  It  thus  behooves  us  to  study  phase  retrieval  in  terms 
of  the  Mellin-transform  power  spectra  as  well  as  spectra  of  the 
Fourier  transform. 
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The  use  of  comparison  filters  in  linear  filter  theory 

Harry  E.  Moses** 

Center  for  Atmospheric  Research,  University  of  Lowell,  Lowell  Massachusetts  01854 
(Received  IS  April  1983;  accepted  for  publication  24  June  1983) 

In  the  present  paper,  it  is  shown  how  the  linear  filter  equation  for  a  given  correlation  coefficient 
can  be  solved  in  terms  of  the  solution  of  the  filter  equation  with  a  different  correlation  coefficient. 
The  second  filter  is  called  a  comparison  filter.  One  obtains  an  integral  equation  for  the  difference 
of  the  two  filters  in  terms  of  the  difference  of  the  two  correlation  functions  and  the  solution  of  the 
comparison  filter.  Thus  if  the  comparison  filter  is  known  and  its  correlation  coefficient  is  close  to 
that  of  the  desired  filter,  one  may  regard  the  comparison  filter  as  being  an  approximation  to  it.  The 
difference  of  the  two  filters  is  then  small  and  perturbation  expansions  or  variational  principles  for 
the  difference  may  be  expected  to  give  better  results  than  if  one  did  not  use  a  comparison  filter. 
The  difference  in  the  solutions  of  the  two  filter  equations  may  also  be  regarded  as  the  change  (or 
error)  in  the  filter  due  to  a  change  (or  error)  in  the  correlation  coefficient.  Our  result  is  obtained  by 
pressing  the  close  analogy  of  the  filter  equation  to  the  Gel’fand-Levitan  equation  of  inverse 
spectral  theory.  Another  result  of  the  use  of  comparison  filters  is  to  show  that  the  filter  equation 
for  the  difference  of  filters  satisfies  a  possibly  useful  grouplike  property. 

PACS  numbers:  02.30.  +  g,  02.30.Rz 


L  INTRODUCTION 

In  Ref.  I,  Kay  and  Moses  treated  the  Gel'fand- Levitan 
equation  of  the  inverse  spectral  theory  problem  from  a  very 
general  point  of  view  and  observed  that  the  Gei'fand-Levi- 
tan  equation  was  a  generalization  of  the  filter  equation  of 
that  time,  namely,  the  Wiener-Hopf  equation.  This  observa¬ 
tion  continues  to  hold  for  more  general  filters,  for  example, 
the  Kalman  filter.  Recently,  one  of  us  (Moses,  Ref.  2)  gave  a 
general  scheme  for  introducing  comparison  potentials  for 
which  the  solution  of  the  corresponding  Oel’fand-Levttan 
equation  is  known.  The  solution  of  any  other  Gel’fand-Levi- 
tan  equation  could  be  expressed  in  terms  of  the  known  solu¬ 
tion  through  the  use  of  an  integral  equation  for  the  difference 
of  the  known  and  sought  for  Gel’fand-Levitan  kernels.  The 
use  of  comparison  potentials  led  to  perturbation  schemes 
and  variational  principles  which,  in  principle  at  least,  led  to 
more  accurate  approximations  for  the  desired  Gel'fand-Lc- 
vt tan  kernel. 

The  purpose  of  the  present  paper  is  to  give  the  analog 
for  the  filter  equation.  It  is  shown  how  the  linear  filter  equa¬ 
tion  for  a  given  correlation  coefficient  can  be  solved  in  terms 
of  the  solution  of  the  filter  equation  with  a  different  correla¬ 
tion  coefficient.  The  second  filter  is  called  a  comparison  fil¬ 
ter.  One  obtains  an  integral  equation  for  the  difference  of  the 
two  filters  in  terms  of  the  difference  of  the  two  correlation 
functions  and  the  solution  of  the  comparison  filter.  Thus  if 
the  comparison  filter  is  known  and  its  correlation  coefficient 
is  doae  to  that  of  the  desired  filter,  one  may  regard  the  com¬ 
parison  filter  as  being  an  approximation  to  it.  The  difference 
of  the  two  filters  is  then  small  and  perturbation  expansions 
or  variational  principles  for  the  difference  may  be  expected 
to  give  better  results  than  if  one  did  not  use  a  comparison 
filter. 

The  difference  may  also  be  regarded  as  the  change  (or 
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error)  in  the  filter  due  to  a  change  |or  error)  in  the  correlat 
coefficient.  Another  result  of  the  use  of  a  comparison  Rite 
to  show  that  the  filter  equation  for  the  difference  of  filter 
satisfies  a  possibly  useful  grouplike  property. 

II.  THE  FILTER  EQUATION  AS  A  GEL’FAND-LEVITAl 
EQUATION 

In  dealing  with  the  filter  equation,  we  shall  use  stands 
notation  as  given,  for  example,  in  Kailath's  monograph  (R 
3).  The  filter  equation  is  then 

h  (M)  =  K  M  -  j*  h  U,r)A  (rjtfr  (/0<j</<// 

In  Eq.  ( 1 ),  h  (/,j)  is  the  filter  matrix  h  (tjt)  —  |  h,t  (r,j)  | ,  and  tj 
matrix  K (te)  =  { Kit(tr r)(  is  related  to  the  signal  correlatit 
matrix  ffr(M)  by 

=  l,&(t  -s)  +  K  MasEXf  ]y'{t  |.  I 

The  Gel’fand-Levitan  equation,  as  treated  in  Ref.  1.  is  ide 
tical  to  the  filter  equation  in  which  A  (/a)  is  the  negative  of  t 
Gel’fand-Levitan  kernel  </  \K  |s),and  the  matrix  AT  (Lr)ist 
driving  kernel  (r  \S1  |.r)  of  Ref.  1.  To  press  the  analogy  ev 
further,  it  is  useful  to  define  the  filter  matrix  as  having  th 
triangularity  property 

A  (1,55  =  0  for  s>t. 

In  the  space  of  observables  which  include  the  signal: 
y{t ),  let  us  define  the  operators  in  terms  of  integral  operate 
with  kernels.  For  example,  if/(r )  is  in  the  space,  we  shall  wr 

A/(r )  £a  Mf[s)ds,  RJ(t )  .  J 

and  so  on.  Let  us  define  the  operator  U  by 
Uf[l  )=/(»)  -A/(l). 

In  particular,  if  y(t )  is  a  signal. 


IMO 
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